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Abstract 

Approximately 10 years ago, the method of renormalization-group symmetries 
entered the field of boundary value problems of classical mathematical physics, stem- 
ming from the concepts of functional self-similarity and of the Bogoliubov renormal- 
ization group treated as a Lie group of continuous transformations. Overwhelmingly 
dominating practical quantum field theory calculations, the renormalization-group 
method formed the basis for the discovery of the asymptotic freedom of strong 
nuclear interactions and underlies the Grand Unification scenario. This paper de- 
scribes the logical framework of a new algorithm based on the modern theory of 
transformation groups and presents the most interesting results of application of 
the method to differential and/or integral equation problems and to problems that 
involve linear functionals of solutions. Examples from nonlinear optics, kinetic the- 
ory, and plasma dynamics are given, where new analytical solutions obtained with 
this algorithm have allowed describing the singularity structure for self-focusing of 
a laser beam in a nonlinear medium, studying generation of harmonics in weakly 
inhomogeneous plasma, and investigating the energy spectra of accelerated ions in 
expanding plasma bunches. 
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1 Introduction 

The paper presents materials illustrating the use and extensions of the concepts of func- 
tional self-similarity and the Bogoliubov renormalization group in boundary value prob- 
lems of mathematical physics. 

The (Lie transformation) group structure discovered by Stueckelberg and Petermann 
in the early 1950s in calculation results in renormalized quantum field theory and the exact 
symmetry of solutions related to this structure were used in 1955 by Bogoliubov and one 
of the present authors to develop a regular method for improving approximate solutions of 
quantum field problems, the renormalization group (RG) method. This method is based 
on the use of the infinitesimal form of the exact group property of a solution to improve 
a perturbative (that is, obtained by means of the perturbation theory) representation of 
this solution. The improvement of the approximation properties of a solution turns out 
to be most efficient in the presence of a singularity, because the correct structure of the 
singularity is then recovered. 

The most spectacular results obtained by the renormalization-group method in quan- 
tum field theory were the discovery of the asymptotic freedom of non-Abelian gauge 
theories (Nobel Prize in 2004), which led to the creation of quantum chromodynamics, 
and sketching the picture of the joint evolution in energy of the three effective interaction 
functions (electromagnetic, weak, and strong) in the Standard Model, which led to the 
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speculative conjecture of a Grand Unification of interactions and the possible instability 
of the proton. 

Apart from this, the quantum field renormalization group provided a foundation (see 
K.G. Wilson's Nobel lecture, 1982) for the construction of an approximate semigroup in 
the investigation of phase transitions in large spin lattices, the so-called Wilson renormal- 
ization group, which is widely used in the analysis of critical phenomena. 

In the present paper, we discuss the most interesting results obtained by the authors 
by extending the RG concepts in quantum field theory to boundary value problems of 
classical mathematical physics. The main achievement here was the development of a 
regular algorithm for finding symmetries of the RG type by means of the modern theory of 
transformation groups. The existence of such an algorithm eliminates the usual deficiency 
of the RG approach beyond the scope of quantum field theory problems: finding the 
group property of solutions requires using special-purpose methods of analysis, usually 
nonstandard, in each particular case. 

It is notable that the algorithm of the construction of renormalization-group symme- 
tries proposed here can be applied to problems involving differential and integral equa- 
tions, as well as linear functional of the solutions. 

We illustrate applications of the algorithm by examples from nonlinear optics, kinetic 
theory, and plasma dynamics, including the problem of propagation and self-focusing of 
a wave beam in a nonlinear medium (Sections 13.21 and 13. 3p . problems of the dynamics of 
a plasma bunch and ion acceleration (Section 14. 2p . and the generation of harmonics in 
laser plasma (Section 13. ip . There, the use of renormalization-group symmetries brought 
about new exact and approximate analytic solutions of nonlinear physics problems, which 
allowed describing the space structure of a self-focusing beam in a nonlinear medium 
in a realistic setting, making significant progress in establishing relations between the 
intensity of the harmonics generated by weakly inhomogeneous laser plasma in a strongly 
nonlinear regime and the parameters of the radiation and the plasma, and finding, for 
the first time, the energy spectrums of accelerated ions in the kinematic description of an 
adiabatic expansion of plasma bunches consisting of several kinds of ions. 

This paper is motivated by our desire to draw theorists' attention to a new and fairly 
general algorithm based on applying the symmetry to an approximate solution for en- 
hancing its approximation power. The use of the group property (the symmetry) of a 
solution underlies both the renormalization group method in quantum field theory and 
its analogue, the new renormalization group symmetry algorithm in mathematical physics. 

The universality of the renormalization-group ideas allows a unified approach to the 
analysis of properties of solutions of various nonlinear problems and gives grounds for 
hopes that this method can be efficiently used in other areas of contemporary physics. 

As is known, this universality is a characteristic feature of another general method 
that represents a solution as a 'path integral' (functional integral) and is widely used in 
quantum mechanics, quantum field theory, the theory of large statistical systems, and 
turbulence theory. 

Classical mathematical physics deals with physical objects described by (ordinary or 
partial) differential equations, which are nonlinear, or integrodifferential in most practi- 
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cally interesting cases. Finding analytic solutions of such equations for arbitrary initial 
and/or boundary conditions is impossible: normally, exact analytic solutions can only be 
found for initial and boundary data of a special form; in other cases, we must content 
ourselves with approximate solutions. The method of constructing a solution of a spe- 
cific boundary value problem (BVP) is usually peculiar to the equations of the particular 
problem under consideration. 

In this paper, we present a method of investigation of analytic solutions based on 
the construction and use of symmetries of a special form of BVP solutions, which we call 
symmetries of the renormalization group kind or renormalization group (RG) symmetries. 
We treat the notion of 'symmetry' in the standard sense of continuous transformation 
groups: this means that a solution of the BVP is transformed into another solution of the 
same BVP by a continuous transformation group acting in the space of all the variables 
determining the solution. The attribute 'renormalization group' points to similarities 
existing between these symmetries and the symmetries in quantum field theory related to 
the operation of renormalization of masses and charges (coupling constants). 

We note that a connection between symmetries and the problem of finding solutions 
of differential equations was first established [1,2] by a Norwegian mathematician, Sophus 
Lie (1842-1899), who showed that most results on the integration of ordinary differential 
equations of various kinds can be obtained by a general method, subsequently called the 
group analysis of differential equations. As one of the main ingredients of the theory of 
continuous groups, the group analysis of differential equations allows classifying differen- 
tial equations using the language of symmetry groups, i.e., it produces a complete hst of 
equations that can be integrated (or such that their order can be reduced) by the group 
method and also suggests a regular procedure for finding these symmetries. Considerable 
progress in this area since the early 1950s has led to new concepts and algorithms, and 
has also extended the range of possible applications of the group analysis (see, e.g., mono- 
graphs [3-9] and handbook [10]), but it has not changed the general aims of the modern 
group analysis to develop regular methods of constructing and classifying solutions of 
nonlinear differential equations on the basis of the symmetries of these equations. 

In problems described by ordinary differential equations, the use of a symmetry group 
yields general and particular solutions. In problems involving partial differential equa- 
tions, which are typical in mathematical physics, knowing a symmetry allows constructing 
particular solutions of a BVP (invariant solutions, which are mapped into themselves by 
the group transformations, and partially invariant solutions), with boundary data not 
known a priori and determined in the construction of a specific solution. Because arbi- 
trary boundary data are not normally invariant under group transformations, the use of 
invariant solutions is generally considered inefficient for the solution of BVPs. 

Arguments underlying the renormalization group method in quantum field theory 
lead to a different conclusion [11]. This method uses the group property of a solution 
(expressed in quantum field theory as a functional equation) for the enhancement of its 
approximation power. 

Although the renormalization group method was originally formulated for quantum 
field problems, we can explain its core idea by an example of a planar problem of radi- 
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ation transfer [12,16]. We assume that the half-space x > is filled with homogeneous 
matter and a stationary stream of particles, characterized by a number a, is falling on 
the boundary a; = of this medium. We consider the evolution of the number of particles 
in the stream as it moves deeper into the medium. Let ai be the number of particles 
in the stream at a distance x = li from the vacuum-matter interface and a2 the number 
of particles at the distance x = I2 = li + \ from the interface. Because the medium is 
homogeneous, the number of particles moving inside at a distance / from the interface is 
uniquely determined by some function A{1, a) of the value of a at the interface and the 
distance /, i.e., ai = A{li,a) and 0^2 = A{l2,a). But the value 0:2 can also be expressed 
as 0:2 = ^(A, ai) in terms of the same function A{X, ai) of two variables, the distance 
A from the imaginary interface x = h and the number of particles ai at this interface. 
Combining the two different definitions of 0:2, we obtain the functional equation 



for A{x,a). The nature of the particles and the properties of the medium are irrelevant 
for this consideration. Of course, solving the transport problem (i.e., an integrodifferential 
kinetic equation), one can find the explicit (exact or approximate) form of the function 
A{x, a) in each particular case, but the exact solution of the problem necessarily satisfies 



A functional equation of form ([T]) occurs naturally in considering a one-parameter 
group G of point transformations in the plane 

Ta : X = f{x,u,a) , u = g{x,u,a), f{x,u,0) = x, g{x,u,0) = u, (2) 

mapping a point P = {x,u) into another point P = {x,u) = Ta{P). We recall that a set 
G of invertible transformations forms a group if these transformations satisfy several 
conditions: a) the set G contains the identity transformation Tq; b) each has an inverse 
transformation T^-i; c) the composition TbTa of two transformations is also an element of 



The last condition in (j3]) can be expressed in terms of the functions / and (7 in ([2]) as two 
functional equations: 

f{x, u,a + b) = f{f{x, u, a),g{x, u, a),b), 
g{x, u,a + b) = g{f{x, u, a),g{x, m, a),b) . 

It is known from the Lie theory that each continuous one-parameter group is fully deter- 
mined by the infinitesimal transformation 





Eqn (P. 



G: 



T„-i(P) = P; ToiP) = P- P =^ niP) = T,Ta{P) = T^+biP) ■ 



(3) 



X = X + a^{x,u) + O(a^) , u = u + ar]{x, u) + 0{c?) 






which is customarily expressed using the infinitesimal operator (or generator) 



'u 



(6) 
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of the group. Finite transformations of a continuous group are uniquely determined by 
the infinitesimal generator by means of the Lie equations, which are the characteristic 
equations for the first-order partial differential equation associated with (P), 

^ = ^(x,m), x\a=o = x, ^ =r]{x,u), u\a=o = u. (7) 
da da 

For the radiation transfer problem under consideration, we have f = x, g = A{a,X), 
a = X, and u = a, and functional equation ([T]) coincides with the second equation in (jlj) 
(the first equation there is an expression of the obvious additive law of the transformation 
of the coordinate x = X), and the group generator is given by 



X = d^ + T]{a)da , r]{a) = dxA{a, X) 



A=0 



In accordance with (^^, to find A at a large distance from the boundary, i.e., for large 
values of the parameter A, we must know the behavior of a = A{a, X) in a thin boundary 
layer, as A — >■ 0, i.e., we must in fact know the derivative of this function at the boundary. 
This information can usually be extracted from an approximate solution provided by 
the perturbation theory. Next, integration of the Lie equations yields formulas for finite 
transformations: 



a 



X 



A. (9) 



Assuming that \E' has the inverse function \1' ^, we find solutions of functional equation 
([T]) in the general form: 



a 



A(a,A) = ^-^(^(a) + A), x = X. (10) 



These constructions are the essence of the renormalization group method. We now present 
two examples of implementing this method. 

We consider a medium absorbing particles in proportion to their number; from the 
perturbation theory, we know the approximate solution 

Apt{a, X) ^ a — uaX , u = const. (11) 

Calculating the coordinate 77(a) = —ua of group generator ([8]) with the help of (fTT!) and 
using it in relations ([9D, we obtain solution (fTOj) . 

^4^3(0;, A) = a exp (— z/A) , (12) 

which is valid in the entire space filled with matter, up to x — > cxd. 

We now assume that the absorption has a nonlinear mechanism, with the absorption 
coefficient proportional to the stream of particles, z^(a) = l3a, where (3 = const. In a thin 
boundary layer, we then have 

Apt{a,X) ^ a- (3a'^X, (13) 
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and the use of f|T3l) in ([9]) yields a solution of ffTOl) in the form of the sum of a geometric 
progression: 



particular, describes the asymptotic behavior of the permeating stream as x — >■ oo. 

The efficiency of the (renormalization) group approach in the above examples shows 
itself in the following fact: using information about the behavior of a solution in a neigh- 
borhood of the vacuum-medium interface, we obtain explicit expressions for the solution 
over the entire interval ^ x ^ oo. We note that if we expand expression f|T^ in a power 
series in the particle number density, i.e., return to the perturbation theory, then in each 
order in n, we obtain expressions increasing in proportion to A", which is a distorted 
representation for the asymptotic form of the solution. The advantage of the renormal- 
ization group method is the recovery of the actual structure of the solution, consistent 
with functional equation ([1]), which is distorted by perturbation theory approximations!!] 

In the case of the planar problem of radiation transfer, the transparency of the renor- 
malization group method is a consequence of taking account of the symmetry properties 
of solutions (i.e., of the functional equation for them) in the actual configuration space. 
The RG transformation of the particle number density in moving deeper into the medium 
is related to a shift in the spatial coordinate. 

Returning to the renormalization group method in its original (quantum field) formu- 
lation [11, 13-15], which is also called the Bogoliubov renormalization groupj^ we note 
that it is based on a functional equation that in the simplest case has the same form as 
([1]) after the substitution x — In t, such that the RG shift transformation of a spatial vari- 
able in transfer theory corresponds to the transferred momenta rescaling in quantum field 
theory; the quantity a is called the invariant (effective) coupling function in this theory. 
In particular, a solution of form (IT^ with A = Inx occurs in quantum field calculations 
in the one-loop approximation. If a more advanced perturbation-theory approximation is 
used, which differs from fll3p by the presence of terms cubic in a, which corresponds to 
the two-loop approximation in quantum field theory, then the RG-improved solution can 
be found from an equation similar to ffTOl) that is unsolvable in elementary functions [17]. 
It is usually solved iteratively, using the one-loop approximation of the RG expression. 

The comparison of the RG-improved solution found in the two-loop approximation 
with the result obtained in the one-loop approximation reveals a characteristic feature 
of the renormalization group method: we can progressively improve the accuracy, which 
is an indication of the stability of the asymptotic behavior of the solution. Similarly, in 
the perturbation theory, we can also take higher-order corrections into account, which 
successively improves the corresponding RG solutions. 

^ We note that formulas and p4)) can also be obtained in transfer theory by other methods, for 
instance, by solving the kinetic equation; however, the method described here, based on the use of group 
differential equations, is the simplest. Moreover, in several important cases, this is the only possible 
method: results obtained with this method are unattainable in other ways. 

■^This term was introduced in [16] to distinguish the Bogoliubov RG from different constructions also 
called renormalization groups in some other areas of physics. They are briefly listed, e.g., in [16, Sec. 3]. 




(14) 



This result, similarly to the previous one, holds in the entire subspace x > and, in 
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Thus, the procedure for the systematic (successive) improvement of the system of 
approximate solutions found in quantum field theory in the perturbation theory with 
respect to a known small parameter is quite similar to the above. This improvement 
of the approximation properties is most significant in the neighborhood of a singularity 
of the solution. In the quantum field context, these are singularities in the infrared 
(see [13,14,18]) and ultraviolet domains. The latter include the most spectacular result 
obtained with the help of the RG method, the discovery of the asymptotic freedom of 
non-Abelian gauge theories [19]. 

The above examples of the use of the renormalization group method for improving and 
refining the approximation properties of perturbative solutions are based on a functional 
equation of the simplest form, with one or two independent and one dependent variable. 
But the number of independent and dependent variables in the problem is often larger 
than this minimal set. 

For example, a version of functional equation ([1]) with x = \nt corresponds to a mass- 
less model with one coupling constant in quantum field theory. We can make this model 
more involved in two ways. First, the number of arguments defining the effective coupling 
a can be increased. For instance, the field model under consideration can contain one or 
several masses (e.g., as in quantum chromodynamics); in that case, the effective coupling 
acquires a dependence on several mass variables with the corresponding transformation 
laws, with the result that the group transformations and the functional equation change 
their form. Second, the number of functional equations can be larger, which corresponds 
to a quantum field model with several coupling constants. This means that we now con- 
sider a group of continuous transformations of independent variables x = {x^, . . . , x"} and 
dependent variables u = {u^, . . . , m™} with infinitesimal operator (jH]) in the space 
and the coordinates of generator are vectors ^ = {^^, . . . , and r] = {rj^, . . . , 77*"}; 
the corresponding contributions to the infinitesimal operator must be understood as the 
result of the contributions of the individual variables. With an increase in the number 
of arguments of the function to be governed by the functional equation and an increase 
in the number of the equations themselves, finding the group property of the solution 
that can be expressed by a functional equation (if we use the original formulation of the 
renormalization-group method [13]) requires a special and often nontrivial analysis in each 
particular case (see, e.g., the discussion in [16,20]); from the algorithmic standpoint, this 
is a deficiency of the RG techniques outside the quantum field theory. 

To overcome this deficiency in extending the RG concepts to problems of mathematical 
physics, another special RG algorithm was developed (see [21,22] and also reviews [23]; [24, 
p. 232]; [25,26]). It has the same aim of finding an improved solution (in comparison with 
the initial approximate solution) as the algorithm of Bogoliubov's RG method, but in 
finding symmetries of a solution of a EVP it uses a scheme of calculations similar to that 
of the modern group analysis. This feature explains the term 'RG symmetry' [23]. 

In this paper, we describe the RG algorithm in mathematical physics and illustrate its 
capabilities by various examples of EVPs. The paper is organized as follows. In Section [2l 
we explain the core ideas of the RG algorithm using the example of the construction of an 
RG symmetry for a solution of a EVP for the Hopf equation. Sections [3] and H] illustrate 



8 



different approaches to the construction of RG symmetries; furthermore, in Section [3l 
we consider several progressively more complicated problems obtained by modifying and 
supplementing the Hopf equations in Section [2J Section H] follows the same logic, but we 
supplement the presentation there with a discussion of nonlocal problems, which are not 
necessarily connected with the ones in Section[3l The scope for possible applications of the 
RG algorithm and a brief list of results obtained with its use are presented in Section [51 

2 The renormalization-group algorithm in mathemat- 
ical physics 

We preface the description of the RG algorithm with the following simple argument. It 
is known that if we treat all the variables (independent or dependent in the standard 
sense) involved in a differential equation and their derivatives (called differential variables 
in group analysis) as independent, then the differential equation can be regarded as an 
algebraic relation for these variables. In the case of one equation, this relation describes 
a 'surface' in the extended space of all the variables involved in the equation (if there are 
several equations, then we speak of a manifold), and each solution of the equation defines 
a 'line' on this surface. The projection onto the {x, u} 'plane' defines a family of curves, 
one of which passes through the 'point' {xq, Uq} corresponding to the boundary condition 
of the BVP in question. 

Transformations of the group G move points on the surface (the manifold) along this 
surface, and therefore the equation preserves its form in the transformed variables and 
each solution of the equation is taken into another solution. A transformation Ta from the 
group G maps a point on the plane {x, u} G M""*""^ into a point {x, u}, and the geometrical 
locus of these points is a continuous curve (a trajectory of the group G) passing through 
{x,u}. The locus of images Ta{{x,u}) is also called the G-orbit of the point {x,u}. In 
the general case, the motion along a group trajectory corresponds to the transition from 
one curve in the family to another, that is, to a 'multiplication' of solutions. 

Returning to the renormalization-group point of view, we consider only the group 
transformations under which points on the curve passing through {a;o, uq} are moved 
along this curve. This means that the solution of the BVP is the RG orbit of the point 
{xo,Uo} (of the boundary manifold in the general case) and is an invariant RG manifold 
(similarly to the invariant charge in quantum field theory [15]). We use the infinitesimal 
version of this property in our construction of the RG symmetry. 

The group property of a solution of a BVP manifests itself as follows: instead of the 
boundary point {xo,Uo} parameterizing the solution, we can take another point in this 
curve related to it by an RG transformation. This 'universality' of the solution of a 
BVP under a change of the way of parameterization is called 'functional self-similarity' 
[27]. To find RG transformations that map a solution of a BVP into a solution of the 
same BVP, we use the fact that a physical problem is formulated in terms of differential 
(integrodifferential) equations whose symmetries can be found by the techniques of group 
analysis. 
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We now illustrate the characteristic features of the algorithm for constructing an RG 
symmetry by an example of a BVP for the Hopf equation [26], which is widely used 
in physics for the description of the initial perturbations at the nonlinear stage of their 
evolution: 

dtv + vdxV = , v{0,x) = eU{x) , (15) 

where U is an invertible function of x and the parameter e defines the 'amplitude' of the 
initial perturbation 'at the boundary' t = 0. For a very small distance t <^ 1/e from the 
boundary, the solution of problem f|T5l) given by the perturbation theory is a segment of 
a power series, 

v = eU- eHUd^U + O (t^) , (16) 

but this form becomes inapplicable for large t. The RG symmetry allows improving the 
perturbative result and recovering the correct behavior of the solution in a neighborhood 
of a singularity (when such a singularity occurs for some values of t). 

In constructing an RG symmetry, the algorithm uses the symmetry group of the 
BVP equations. The boundary data defining a particular solution are involved in RG 
transformations by extending the space of the variables on which the group acts. In the 
case of BVP ffT5]) . this space involves three independent variables, x = {t,x,e}. It is 
convenient to write differential equation f[T^ for the function u = v/e introduced such 
that the 'amplitude' e is carried over from the boundary condition to the differential 
equation: 

dtu + eudxU = , ^(0, x) = U (x) . (17) 

The general element of the transformation group G for Eqn f|T7j) (for the basic manifold 
in the general case) can be found by means of the standard Lie techniques (see, e.g., [4]); 
it is given by a combination of four infinitesimal operators, 

(18) 

X3 = ijj^ {xdx + udu) , X4 = 7/;^ {ed^ + xd^) , 

where ip'^ {i = 2,3,4) are arbitrary functions of e, u and x — eut and ip'^ is an arbitrary 
function of all the group variables {t,x,e,u}. We now use the RG invariance condition 
for a particular solution of BVP f|T7|) defined by the relation 

S = u-W{t,x,e) = (19) 

with the function W that is unknown at this point; in other words, we check that the RG 
transformation maps the solution of the BVP into the same solution. In the infinitesimal 
form, this condition can be written as 

X5i ^, = ip%W - xdxW) - ip^dxW - %l)\ed,W + xd^W) = , (20) 

where \[s] means that the result of the action of the operator is taken on the manifold 
defined by the equation S = and all its differential consequences. The term containing 



10 



is absent in fl20l) because it is proportional to dtW + eWd^W, which vanishes identically 
on solutions of Eqn f[T7|) . Condition fl2U]) holds for all t, and for t — > in particular, when 
W is replaced by the approximate solution 

W = U- etUd.^U + O (t^) (21) 

obtained in the framework of perturbation theory (fT6|l . In this limit, Eqn (120|) yields a 
relation for the functions (i = 2, 3, 4), which extends in the obvious fashion to t 7^ 0: 

V^' = + V^') + (w/^xf^)^' , X = ^ - , (22) 

where the derivative d^U must be expressed in terms of x or w in accordance with the 
boundary conditions. Using (l22l) in ( JTSll . we arrive at a group of a smaller dimension with 
the infinitesimal operators 

i? = Vi?,, i?i = 7/-^ (9t + ewa,) , 

(23) 

R2 = u^j^ [(et + lld^U) 9, + 9„] , i?3 = e^' (tM4 + d,) . 

The above procedure reducing (ITSl) to (1251) is t/ie restriction of group f ITS]) on a particular 
solution^ and the set of operators i?j in (l23l) describes the required RG symmetry. We 
obtain the solution of the BVP with the use of the corresponding Lie equations (similar 
to ([7])) for any generator in (!23|) . Without loss of generality, we can take the generator 
i?3 with eip^ = 1 to obtain the finite RG transformations 

x' = X + atu , e' = e + a , t' = t , u' = u , (24) 

where a is the group parameter, t and u are invariants, and the transformations of e and 
X are translations, which in addition depend on t and u for the x variable. For e = 0, in 
view of f|T7|) . the variables x and u are related hj x = H{u), where H{u) is the function 
inverse to U{x). Eliminating a, t, and u from flM|) and dropping the dashes in our notation 
for the variables, we obtain the required solution of BVP f[T7|) in implicit form [similar to 
the implicit form of the solution of functional equation flTU]) ]: 

x-etu = H{u) . (25) 

In effect, this is the improved perturbation theory solution (fT6|) . which can be used not 
only for small t -C 1/e (of course, under the condition that fl25|) defines u uniquely). 
Depending on H{u), this solution either indicates the correct asymptotic behavior as 
t ^ 00 or gives the correct description of the solution in the neighborhood of finite values 
t — ^ tging. One example of the first option is the solution of the BVP for the linear 
function U{x) = x. This yields the expression v = ex{l + et)~^, which remains finite 
as t — >■ 00, similarly to the solution of f|T^ . For the second option, we can select, for 
instance, a sine wave U{x) = — sin x at the boundary. Then solution (|25|) describes the 
well-known distortion of the initial profile of a sine wave, transforming it into a saw-tooth 
shape [28, Ch.6, §1], with a singularity forming at a finite distance tging = 1/e from the 
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boundary. We note that for finding solution fl25l) of the BVP, we use only the known 
symmetry of the solution and the corresponding perturbation theory (PT). 

The above example of the construction of RG symmetries illustrates the general al- 
gorithm, whose detailed description in relation to EVPs for differential equations can be 
found, e.g., in reviews [23,24], and whose generalization to nonlocal problems is presented 
in [25,26]. We can schematically express the implementation of the RG algorithm as a 
sequence of four steps (see the figure): 



(I) constructing the 

basic manifold TZA4; 

(II) finding a symmetry group G 
admitted by TZJ^; 

(III) restricting the symmetry 
group G on a particular solu- 
tion of the BVP and finding 
the RG symmetry (RGS); 

(IV) finding an analytic solution 
corresponding to the RG sym- 
metry. 



- PT solution 



Basic model 



Basic manifold 1 



Generators of the group G 



III 



— > Generators of the RG symmetry 



IV 



RG-invariant solution 



Figure. The scheme of the RG algorithm 



A characteristic feature of the procedure of constructing the RG symmetry is the mul- 
tivariance of step (I) , whose aim is to have the parameters participating in the equations 
and the boundary conditions of the problem and determining the solution somehow in- 
volved in transformations. The choice of a concrete realization of the first step is most 
usually governed by the form of the basic equations and the corresponding boundary con- 
ditions on the one hand and by the form of the approximate PT solution on the other. 
This multivariance, which is a feature of step (I) alone, is aimed at covering a possibly 
broader spectrum of problems to be investigated by the method. The subsequent steps 
are carried out in the framework of well-developed group methods. 

This multivariance is also seen in the above simple example of a BVP for the Hopf 
equation. Underlying our construction of the RG symmetry for BVP (1171) was the most 
obvious option: constructing the RG symmetry from the point symmetry group of the 
Hopf equation in the space extended by incorporating the parameter e into the set of 
independent variables. This way of constructing the basic manifold is not the only 
possible one. 

We could also construct the RG symmetry for BVP flTTl) using an additional differ- 
ential constraint compatible with the boundary conditions and the basic equationsjf] For 

^Here, we do not detail the construction of such a differential constraint. As an example, we note the 
use of the invariance condition for the basic equation under so-called higher (or Lie-Backlund) symmetries 
rather than point symmetries. In contrast to the coordinates of an infinitesimal generator of a point 
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instance, if the initial conditions in f|T7|l are given by the hnear function U {x) = x, then 
we can take the differential constraint d^xU = 0. Next, we calculate the RG symmetry of 
BVP f[T7l) taking the basic manifold TZAi to be the system obtained by combining this 
constraint and the Hopf equation. The admissible group G for the manifold T^A^ in so 
doing is different from (fTSll . but the form of solution (125|) is the same. Other examples of 
the implementation of step (I) of the algorithm can be found in [23]. 

3 Renormalization-group symmetries in local prob- 
lems of mathematical physics 

The example of the construction of the RG symmetry for Hopf equation f[T7|) demonstrates 
that a particular form of the realization of the general scheme of the RG algorithm de- 
pends on the form of the equations in the BVP, as well as on the way the boundary 
data are specified. Since the construction of the RG symmetry proceeds by restricting 
the symmetry group G of the basic manifold [step (III)], the RG usually has a smaller 
dimension than G. For instance, in the case of BVP (1171) . the symmetry group G is de- 
fined by the four generators Xj, and the RG is defined by the three generators Ri. It is 
obvious that for the construction of the RG symmetry, it is desirable to have a maximal 
group G. However, the more complicated the basic equations are, the narrower the ad- 
missible transformation group typically is. For instance, if the term lyu^x accounting for 
dissipation is added to the Hopf equation, then after the change of variables = w, we 
obtain the modified Burgers equation. For this equation, the admissible symmetry group 
is infinite-dimensional, but it is now characterized by a single arbitrary function instead 
of four functions for the Hopf equation, and after the reduction procedure, we obtain a 
finite-dimensional (8-dimensional) RG [29]. 

It is also possible that the RG symmetry cannot be constructed using a point symmetry 
group for the basic manifold alone because restricting on a particular solution yields a 
zero-dimensional group. In this case, we must either modify (and simplify) the system 
of equations used for the description of the physical process or use other symmetries in 
addition to Lie symmetries for constructing the RG. 

We now demonstrate various approaches to the construction of the RG symmetry for 
the BVP obtained by complicating the problem in f[T7|) . which was our example of the 
construction of the RG symmetry in Section [21 

symmetry group, the coordinates of a generator of a higher-symmetry group in addition to independent 
and dependent variables, also depend on higher derivatives. Expressing the invariance condition under 
a group of higher symmetries in the infinitesimal form, we obtain the required differential constraint 
(see [23] for the details). 
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3.1 Renormalization-group symmetry in nonlinear plasma the- 
ory 

We consider the following problem, which was historically the first example of a successful 
application of the RG algorithm. This is the interaction of p-polarized electromagnetic 
radiation with a frequency uj and a 'moderate' (by today's standards) intensity, with 
inhomogeneous plasma [21]. This interaction is described by a system of 2-dimensional 
nonstationary differential equations (the equations of the coUisionless hydrodynamics of 
electron plasma with a self-consistent electromagnetic field) for six functions: the compo- 
nents Bz and E^, Ey of the magnetic and the electric fields, two components Vx-, Vy of the 
velocity of the electrons, and their density n; these functions depend on three variables: 
the coordinates x and y and time t. Our aim is to obtain an approximate analytic solu- 
tion of this system of equations in an arbitrary order of nonlinearity, without confining 
ourselves to the perturbation-theory framework. 

For an arbitrary ion density function n'(x), the basic system of equations admits 
only a finite-dimensional point transformation group, the group of translations along the 
t and y axes. If the ion density is a constant, if = N = const, then we also have the 
group of translations along the x axis and the group of simultaneous rotations in the three 
planes defined by the coordinates {x, y} and the corresponding x and y components of the 
velocity of the electrons and of the electric field. Thus, regarding the original equations 
as the manifold we obtain a fairly narrow admissible group, which does not allow 

finding the required RG symmetry. 

To construct a manifold TZAi allowing a wider point transformation group, we use the 
fact that the leading contribution to nonlinear effects of the interaction of the electromag- 
netic wave with the inhomogeneous (in x) plasma under consideration here comes from the 
components of the electric field and the velocity of the electrons that are directed along 
the density gradient. Furthermore, due to the natural smallness parameters (the smooth 
inhomogeneity of the ion density along the x axis and the small angle of incidence 6 of 
the laser beam on the plasma), the dependence of these components on the y coordinate, 
which is transverse to the density gradient, is smoother than their dependence on x in 
the neighborhood of the plasma resonance. Hence, in the construction of TZAi, in the full 
system of 6 original equations, we can single out a simpler system of two one-dimensional 
nonlinear partial differential equations for the x components E^ of the electric field and 
Vx of the velocity of the electrons in the neighborhood of the plasma resonance: 

LudrV + avdxV — p = , ojdrP + avdxP + uj\v = , t = ujt — (ojyj c) sin Q . (26) 

Here, v and p are the respective quantities Vx and Ex normalized by the parameter a, 
the parameter a oc ^Jq is determined by the radiation flux q on the plasma and the linear 
transformation coefficient, lol{x) is the plasma frequency (for the fixed ion density), and 
c is the speed of light. 

The infinite-dimensional point transformation group in the space of 5 variables {r, x, a, 
admitted by fl26|) is defined by an infinitesimal operator, which is a sum of three 
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operators: 



1 1 

tr « « (27) 

X3 = —{ada — vdy — pdp) , y = udr + avdx + pdy — uj\vdp . 

(X 

Each of these three operators involves an arbitrary function /ij of the group variables, 
where ^2 and /is satisfy the differential constraints 

Y\ii2) + Y{ul^i2) = 0, F(/i3) = 0. (28) 

To specialize the functions /xi, /i2, and /is entering the coordinates of the operator X, we 
use the procedure of restriction of the point transformation group fl271) on an approximate 
(in powers of a) solution of the BVP. We can construct this solution such that the zeroth 
approximation to the functions v and p is found by solving the linearized system of the 
original six equations endowed with the corresponding boundary conditions (an electro- 
magnetic wave falling on plasma from the vacuum) and by the selected density profile 
n*(x) in the plasma resonance region; corrections to this solution that are proportional 
to a arise after the linearization of system (l26l) . The verification of the RG-invariance 
conditions [similar to (1201) ] for this approximate particular solution determines our choice 
of the functions /ii = 0, 1^2 = —p/^'^i and /is = 1 and yields the required RG-symmetry 
operator (where the first relation in fl28|l holds with the substitution uj\ cj^): 

i? = X2 + Xs = -{pluo^)dx + da. (29) 

The quantities r, v and p are invariants of the RG transformations with infinitesimal 
operator fl29|) . and the transformation of the x variable defined by the solution of the Lie 
equation for fl29|) exhibits a linear dependence on the parameter a: 

X = ri — [p/ uj'^)a. (30) 

The group composition law for x can be easily deduced from the functional equation of 
form ([T]) with the substitutions A — > a;, / — > 02, A ^ ai, and a — > r/. We note that in 
contrast to the transfer theory problem, the group parameter here is not an independent 
variable involved into the equation but the parameter a imported into the equation from 
the boundary conditions. 

The solution of Eqn fl26|) constructed with the help of (!29l) is given by 

ap/t^^A = -£:(/i sin r + /2 cost) , e = {q/qof''^, ^^^^ 
av/uA = e (/i cos r — /2 sin r) , x = rj + e {fi sin r + f2 cos r) , 

where the parameter e oc a oc ^yq, which depends on the flux go of the plasma wave 
breaking at the critical point, does not exceed 1, and the functions fi,2{v) are determined 
by the well-understood linear structure of the field, whose explicit form can be various. 
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depending on the density profile and the thermal motion of the electrons in the plasma. 
In cold plasma with a linear density profile, we have 



/i=(l + r/^)"\ h = ri{l+vT'- (32) 

When a weak thermal motion of the electrons is taken into account, relations (!32l) must 
be modified: 

oo oo 

h = j cos(r/e + eV3) , f2 = j sin(r/e + ^'/S) • (33) 







Solution fl3Tl) is an exact solution of Eqns fl26|) for ujl = uj. The x and t] variables in 
relations fl3T|) - fl33|) . in view of the normalization by the width of the plasma resonance A, 
are dimensionless quantities. The equations for the remaining four normalized quantities 
(the electric field Ey, the magnetic field i?^, the y-component Vy of the velocity of the 
electrons, and the density n) are given by 

dr^Ey = -(cuA/c) sm6drE:^ , LodrVy = Ey , 

d,B, = {Vjc)dr,Ey-{Vy/c)d,E,, n^uj\l/d,x). ^ ' 

Integration of Eqns flM|) is elementary. Formulas fl3T|) and fl34l) present the required 
solution of the BVP. Discarding strongly nonlinear effects, we can use f l3T|) and flM|) to 
obtain results from the theory of generation of arbitrary-order harmonics in cold [30] and 
hot [31] inhomogeneous plasma (if we respectively use formulas fl5^ and (1551) for /i,2)- 
Taking strong nonlinearities (the influence of higher harmonics on the lower ones) into 
account significantly changes the dependence of the coefficient of the transformation into 
harmonics emitted by the plasma on the density of the electromagnetic radiation flux 
falling on the plasma [21,32] and the temperature of the plasma [33,34]. 

Result (15T]) . of solving the BVP for the six original equations takes both the 
boundary condition and the strongest nonlinearity into account, and is exact in the same 
measure in which the group symmetry of Eqns fl26|) reflects the symmetry of the full 
system of six original equations under the above assumptions. The approximate nature of 
the group with infinitesimal operator fl2Ul) so obtained relative to the group fl27|) inducing 
it is determined by the inhomogeneity of the plasma (we recall that in the derivation of 
operator (l27Il . we imposed no assumptions on Eqns (126|) concerning the inhomogeneity 
pattern of the plasma density). This is similar to the situation in quantum field theory: 
the exact group property of a solution is used for a progressive improvement of the system 
of its approximation characteristics, where the next approximation improves the previous 
one without destroying it. From the standpoint of the RG symmetry, this means that 
operator fl29|) can be refined by accounting for the small parameters of the problem used 
in passing to Eqns fl26l) . We say in this case that the symmetry of system fl26l) is inherited 
by a more general system of equations. An example of an RG symmetry for fl26l) with 
the corrections due to the inhomogeneity of the plasma taken into account is presented 
in [21]. 
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3.2 Renormalizat ion-group symmetries in problems of gaseous 
and quasi- Chap lygin media 

The situation where the existence of an infinite-dimensional point transformation group 
ensures the construction of an RG symmetry, as in the examples of BVPs for the Hopf 
equations and Eqns fl26|) . is not universal. Below, we present an example of a EVP in 
which the symmetry group for the original manifold (the system of differential equations) 
is infinite dimensional, but the construction of the RG symmetries requires using higher 
symmetries (which are also called Lie-Backlund symmetries [5]) instead of a point trans- 
formation group. 

We consider the BVP for a system of two nonlinear first-order partial differential 
equations for functions v and n > 0: 



with constant a and a nonlinearity function (f depending only on n. Depending on the 
sign of aip{n), these equations are of either the hyperbolic ( a({>{n) < 0) or the elliptic 
{aip{n) > 0) type. In the first case, fl35l) corresponds to the standard equations of gas 
dynamics for one-dimensional planar isoentropic motion of gas with the density n and 
velocity v. The second case relates to equations of quasi- Chaplygin mediaQ 

Because Eqns (135!) are linear in the hodograph variables t = nt and x = 2; — f t. 



there exists an infinitesimal operator of an infinite-dimensional subgroup Xoo = C,^dr + 
^"^d^, whose coordinates and are defined by arbitrary solutions of the partial differ- 
ential equations transformed into (l36l) by the substitution t — x This means 
that we can formally construct an RG symmetry by restricting an infinite-dimensional 
point transformation group, but this requires knowing the solution of (1361) for arbitrary 
boundary data; in fact, such a procedure is equivalent to solving the original BVP. 

Another approach to the construction of RG symmetries for the problem under con- 
sideration is the use of a higher (Lie-Backlund [37]) symmetry group. By contrast to 
point transformation groups with generators of form ([6]), Lie-Backlund symmetries are 
characterized by an infinitesimal operator with the coordinates depending on indepen- 
dent variables x and differential variables u and the derivatives = {d^iu"} = {u"}, 
U(2) = {dxixju"} = {ufj}, . . ., where a = 1, . . . ,m; i, j = 1, . . . ,n. The relation between 
these variables can be expressed by means of the total differentiation operators Di, as the 
following system of equalities: 



ut = A (w") , = D, «) = D,D, (m") , . . . , a = + <du^ + M° 5,. + . . . . (37) 



^The term 'quasi-Chaplygin media' is used in the discussion of nonlinear phenomena developing in 
accordance with the mathematical scenario for the Chaplygin gas, i.e., the gas with a negative adiabatic 
exponent. At first glance, such a model looks like the standard model of gas dynamics, but it corresponds 
to the negative first derivative of the 'pressure' with respect to the 'density.' A characteristic feature of 
quasi-Chaplygin media is a universal mathematical form of various nonlinear effects accompanying the 
development of an instability. 



dtv + vdxV = aip{n)dxn, dtn + vd^n + ndxV = , 
v{0, x) = aW{x) , n{0, x) = N{x) , 



(35) 



d^T - %lj{n) dnX = , d^x + dnT = , = n/aip, 



(36) 
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The Lie-Backlund group theory allows restricting to only canonical operators, which leave 
all the independent variables invariant. This is important, for instance, in the analysis of 
symmetries of integrodifferential equations and in the construction of the RG symmetries 
in problems involving nonlocal equations. For BVP (!35l) under consideration here, it 
is convenient to write the Lie-Backlund group generators for the equations expressed in 
hodograph variables (136|1 : 

X = J2c^X, = Y^ C^ [hdr + g^^y,) . (38) 
i i 

The coordinates fi and gi of generators (!38|) . which are linear functions of the differential 
variables [38], are connected by a system of recursion relations 

where the entries of the matrix recursion operators Lk are linear functions of the operator 
Dn of total differentiation with respect to n. The number of operators Lk depends on the 
form of the nonlinearity function (p{n); in the most typical case ip{n) = n{n + bY, where 
b, i = const, there are three operators: k = 1,2, 3. The action of the recursion operators 
on the coordinates /i = r and gi = x of the physically 'obvious' dilation operator in 
the space of the hodograph variables t = nt and x = 2; — f t yields three operators with 
coordinates fi and gi {i = 2,3,4) linearly depending on the derivatives r„ = (9„r and 
Xn = dnX'-i they are therefore equivalent to infinitesimal operators of the point group. The 
action of the recursion operators Li, . . . , L3 on the first-order symmetries fi, gi {i = 2, 3, 4) 
generates five operators, whose coordinates in the hodograph variables are linear functions 
of these variables and their second-order derivatives. These are Lie-Backlund symmetries 
of the second order. Repeating this procedure several times, we obtain 2s -|- 1 symmetries 
of a fixed order s [38] . 

The infinite system of operators fl38l) (obtained at step (II) of the RG algorithm) 
for Eqns fl5^ (treated as the TZAi manifold) enables constructing the operators of RG 
symmetries and finding the corresponding RG-invariant solutions. The reduction of the 
Lie-Backlund group (step (III) of the RG algorithm) reduces to the verification of the 
invariance conditions / = and g = (similar to (120|) . but generalized to the case of 
Lie-Backlund symmetries) for a concrete solution of the BVP, where the functions / and g 
are arbitrary linear combinations of some coordinates /, and gi of the canonical operators 
of the group and are chosen so as to satisfy the prescribed boundary conditions at t = 0. 
As examples, we give the values of the coordinates of two second-order operators of the 
Lie - Backlund RG symmetry. 
Example 1. 

f = 2n{l - n)Tnn - UTn - 2nv{Xn + nXnn) + nv'^Tnn/2 , ^^^^ 
g = 2n{l + {2- 3n)Xn + V {2nTnn + Tn) + {v'^/2) {uXnn + Xn) ■ 

Example 2. 

f = -n^ lnrar„„ - (ra/2)r„ + r/2 virC'Xnn + (3/2)n\„) , ^^^^ 
g = -n^lnnx„„ + (n/2)(l-h41nn)x„ + x/2-ht;(nr„„-hr„/2) . 
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The operator R with coordinates fHOj) corresponds to the solution of the BVP for Eqns fl35l) 
with a = 1, ^p{n) = 1 for W{x) = and N{x) = cosh~^(a;), and the operator R with 
coordinates fHT]) corresponds to the solution of the BVP for Eqns (1551) with a = —1, 
= 1/n for = and N{x) = exp(— x^). To solve the BVP using RG symmetries 

fHOl) and (HTj) . we must add the invariance condition f = g = to the basic 7?.A1 and 
solve the resulting system of equations (step (IV) of the RG algorithm). 

For RG symmetry fl40p . a solution exists on a finite interval ^ t ^ tsing, until a 
singularity occurs on the axis x = at t = tsmg = 1/2, when drcv{tsing,0) — > oo and the 
value of n remains finite, n{tsing, 0) = 2: 

V = — 2nttanh(x — vt) , n^t^ = ncosh^(x — vt) — 1 . (42) 

From the physical standpoint, solution which was previously obtained in [39], de- 
scribes the evolution of a planar light beam in a medium with a cubic nonlinearity (a 
quasi- Chaplygin medium) for the boundary condition A^(x) = cosh" (x). The quantities 
n and v define the intensity and the eikonal derivative of the beam. 

For RG symmetry (HTj) . the solution describes a monotonic evolution (decrease) with 
time t of the density n ^ 0, while the particle velocity continues to be linearly dependent 
on the coordinate: 

V = xV2qe-'^'/^ , n = e'^'/^ exp (-x^e-^') , t = (v^/2)erfi (q/V2^ . (43) 

Solution fj43|) . which was discussed in [40], describes an expanding plasma layer with the 
initial density distribution N{x) = exp(— x^). 

These two examples demonstrate that by using the Lie-Backlund RG symmetry, we 
achieve the same goals as with point RG symmetries: we can give an adequate description 
of the structure of the solution in the presence of a singularity or can find its asymptotic 
behavior. Although we found RG symmetries (HTl) and (l42l) for the already known so- 
lutions, the RG approach reveals the group structure of these solutions. Previously, to 
obtain these results, the authors imposed some a priori assumptions about the structure 
of the solution. In [41], the reader can find an example of the solution of a BVP with the 
help of Lie-Backlund RG symmetries for 0351) with the initial condition of a more complex 
type, not representable in terms of elementary functions, when the intensity distribution 
of the light beam at the boundary has the form of a smoothed step function. 

3.3 Approximate renormalizat ion-group symmetries in problems 
of quasi- Chaplygin media 

Constructing an RG symmetry on the basis of higher symmetries is justified if the equa- 
tions defining an RG-invariant solution can be investigated analytically. The complexity 
of differential equations usually increases with their order. Hence, the use of higher-order 
Lie-Backlund symmetries in the invariance conditions of the RG symmetry can often limit 
the potential for applications of such symmetries in the case of arbitrary boundary data. 
On the other hand, a restriction on the order of the allowed symmetries narrows the va- 
riety of approaches to the construction of RG symmetries for arbitrary boundary data. 
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For instance, for BVP fl35|) . the symmetry group of the original manifold fl36|) allows only 
2s + 1 symmetries of a fixed order s, which for small s can be insufficient for the construc- 
tion of the RG symmetry for arbitrary N{x). For the extension of the symmetry group 
of the original manifold, we must use the technique of approximate symmetries [35]. 

The central idea here is the use of natural smallness parameters (which we distinguish 
from the parameter with respect to which we construct the PT approximation to be used 
in RG transformations), which are involved in some form in most physical problems and 
which enter the equations as coefficients. For instance, the coefficient a of the nonlinearity 
function (p{n) in fl35|) is such a parameter. The presence of natural small parameters 
allows expressing the required symmetry as a power series in these parameters and taking 
finitely many terms of this series. If we discard the small parameters altogether, then 
the equations defining the are simpler than the original equations and allow a wider 
transformation group, and hence there can be more approaches to the construction of 
the RG symmetry for arbitrary boundary data. An essential point here is the possibility 
to successively account for corrections to the obtained RG symmetry for the system of 
differential equations of the simplified manifold: when this can be done, we say that we 
have constructed a symmetry inherited in a given order in the small parameter. 

We demonstrate how approximations to the RG symmetry for BVP (153]) can be con- 
structed for small a <^ 1. Setting w = v/a, we write system of equations as 

dy,T - {n/ip{n))dnX = 0, dy,x + adnT = . (44) 

As a — > 0, dropping the second term in the second equation yields a simpler subsystem of 
differential equations, which is an approximation to the original manifold TZAi. By con- 
trast to the symmetries of Eqns fl36|) . which allow only a finite-dimensional Lie-Backlund 
symmetry group of a given order, Eqns (jHj) with a = have an infinite-dimensional 
symmetry, which is consistent with the perturbation theory for the BVP with arbitrary 
boundary data. Hence, we seek RG symmetries by combining symmetries of the 'zeroth' 
approximation to the equations (i.e., of Eqns (jUj) with a = 0) and corrections to them in 
powers of a. We represent the coordinates / and g of the canonical operator of the group 
for fH4|) as a power series in a: 

oo oo 

X = fd, + gd^, f = Y^a'f- g = Y^a'g\ (45) 

i=0 1=0 

Using the techniques of modern group analysis [10] that generalize Lie's algorithm to 
higher symmetries, we obtain a system of recursive relations for the /* and g'^: 



gi = G' + ^l- J aw {Zg'-' - Yf-'} 



(46) 
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where 



oo 



oo 



Y 



w = v/a , 



s=0 s=0 



(47) 



dn- 



d' 



T 



Xs = 




p J dn'P 



Xs—p+l ! 



F^{n,Xs,T^s) and Xs, ^s) are arbitrary functions; the integrands are expressed in 

terms of f^, Xs, and w. It is an immediate consequence of fH^ that for small a, the 
symmetry of the equations of the 'zeroth' approximation is inherited by system (jH]) in 
any finite order in a: the corrections do not destroy the symmetry f^,g^ of the 'zeroth' 
approximation. The form of the inherited symmetry (i.e., the expressions for / and g) 
is fully determined by relations ( l46ll : it can be a point symmetry or a Lie-Backlund 
symmetry. 

Because the dependence of the functions and G" on their arguments can be ar- 
bitrary, we can construct RG symmetries for the EVP with arbitrary boundary data: 
the restriction of approximate group fH6|) on solutions of the EVP (step (III) of the RG 
algorithm) is performed, similarly to the case of the exact Lie-Eackhand RG symmetry 
in subsection 13. 2^ by verifying the condition f = g = for a concrete solution of the 
EVP. Here, we choose the functions and so as to satisfy the prescribed boundary 
conditions at t = 0. In particular, for EVP (1351) with W{x) = 0, we can set and 
(i ^ 1) equal to zero and can ensure the boundary conditions by selecting F° and G^. 

We now present two examples of RG symmetries constructed with the use of relations 
(1461) . The first example is related to the EVP for Eqns (153]) with ip{n) = 1 for W{x) = 
and N{x) = cosh~^(x). With these conditions, we can take the following functions /° 
and g^: 



f = 2n{l - n)T2 - nn - 2nw{xi + nx2) , ^° = 2n(l - n)x2 + (2 - 3n)xi • (48) 
Substituting /° and g^ in fH6|) . we find the next terms of series (H5|) . the functions and 



and the substitution of and g^ in (l46l) gives zero values for all /* and g'^ with 2^2. This 
means that the RG symmetry can be expressed in this case by binomials / = /° + a/^, 
g = g^ + ag^i that is, infinite series fHS]) terminate and turn into finite sums, and the 
binomial expressions for the RG symmetry are exact and hold for arbitrary values of a. 
In particular, setting a = 1, we arrive at relations fHOj) . 

For arbitrary boundary data, the infinite series in fH5l) do not automatically terminate, 
and taking only finitely many terms of the series means that the RG symmetry constructed 
with the use of (1451) and fHB]) is approximate in the sense described in [42]. The second 
example corresponds to an approximate RG symmetry for EVP (!35|) with ^p{n) = 1 for 
W{x) = and N{x) = exp(-x^): 




/I = nw^T2/2 , g' = w {Inr^ + n) + {w" 12) (nxa + Xi) , 



(49) 



/ = 1 + 2nxXn + a (-2rr„ + rV^) , g = -2a (rx„ + X^n) ■ 



(50) 
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Here, we omit all the contributions to / and g proportional to the higher powers a* with 
i ^ 2. 

The above constructions of RG symmetries can be easily generalized to the case where 
the group transformations involve the parameter a in addition to the 'natural' variables 
of the problem. In this case, the set of possible RG symmetries is usually larger. For 
example, we note the approximate RG symmetry for the same EVP as in the second 
example, but, in contrast to (ISO]) , containing derivatives with respect to the parameter a: 

f = 2n(rx„ + TnX) + '^(^X'^a , fl- = 1 + 2nxXn + 2a {xXa - rT„) . (51) 

Unlike exact RG symmetries, which allow finding an exact solution of the EVP for any 
RG generator chosen, approximate symmetries yield a solution of the EVP depending 
essentially on the form of the RG symmetry operator, as can be seen, for instance, from 
the use of generators (1501) and (1^ (see [41]). 

The use of several approximate analytic solutions or the comparison of the solution 
obtained on the basis of the exact RG symmetry (and used as a test) with the solution 
obtained on the basis of an approximate RG-symmetry allows evaluating the accuracy of 
the corresponding approximate RG- invariant solution [43]. 

For finding approximate RG symmetries in a physical problem, we can use not one 
but several small parameters. This is the case, for instance, in the construction of the RG 
symmetry for a EVP for the system of equations of a light beam in a nonlinear medium, 
which can be regarded as a natural generalization of (l35l) : 



vt + vv^ - aLp{n)n^ - {38^ ((x " / \/n) dr, {x^dr, {Vn))) = , ^^^^ 
rit + nvx + vux + pinv /x) = , f (0, x) = V{x) , n(0, x) = N{x) . 

The parameters a and f3 determine the contribution of the nonlinear refraction and 
diffraction processes; z/ = for a planar (2-dimensional) wave beam and u = 1 for a 
3-dimensional (axially symmetric) wave beam. 

The construction of an RG symmetry for EVP (I52p proceeds in accordance with a 
scheme similar to the one used before: the coordinates / and g of the canonical operator 
for the manifold TZAi defined by Eqns fl52|) can be represented as double power series in 
the nonlinearity parameter a and the diffraction parameter (3: 

oo oo 

X = fd.+gdn, /=J]aW^ <7=^aW^ (53) 

i,j=0 i,j=0 

The standard techniques of group analysis are used for the calculation of the coefficients 
anci Restricting ourself to finitely many terms of series (l53l) . we arrive in the 

general case at an approximate symmetry, which after the restriction procedure gives 
the required RG symmetry. As an example [44], we present explicit expressions for the 
coordinates / and g of the infinitesimal RG symmetry operator for EVP (1521) in the case 
of a coUimated cylindrical (z/ = 1) beam in a medium with cubic nonlinearity {ip = 1): 

f = Dx!^S-(^an + -^D. (x^xv^)^ | , 9=^0. {{xn) [v - tS^]} , (54) 
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where the function S depends on x = x — vt: 



Six) = aNix) + ^J^ d, [xd^ (x/iVOO)) ■ (55) 

The canonical RG operator with coordinates (l54ll is equivalent to the following operator 
of a point RG symmetry: 



R= [1 + t'^S^^] dt + S^d^ + [tS^ + vfS^^] - nt 



X J X ^ 



dn , (56) 



which allows easily finding finite group transformations (step (IV) of the RG algorithm) 
relating the values of n (the beam intensity) and v (the eikonal derivative) for t > to 
similar quantities at the boundary t = of the nonlinear medium, i.e., constructing the 
required solution of BVP ([52]) [44]. 

In the derivation of fl55l) and fl56l) . we considered contributions of the form /° = /(*^''^) 
and = in fl53|) . that is, contributions independent of a and P, and also con- 

tributions linear in these parameters = af^^'^^ + Pf^^'^^ and = ag^^'^^ + Pg^^'^\ 
Dropping the terms proportional to O (a^,/3^,a/5) means that symmetry fl5B]) is approx- 
imate with respect to these parameters. Of course, similarly to BVPs for equations of 
quasi- Chaplygin media, there exist distributions N{x) for which series (!53l) terminate and 
become finite sums. Such a situation corresponds to the exact RG symmetry in (15^ 
rather than an approximate one and to an exact solution of the BVP for arbitrary values 
of the parameters a and /3. In particular, symmetry fl53p is exact when S{x) is a binomial: 
S{x) = So + S2X^/2. This form of S{x) corresponds to a particular dependence on the x 
variable of the beam intensity N at the interface. For instance, for S2 = and sq > 0, 
Eqn fl56l) yields a solution of the BVP that describes a Townes' self-channeling beam [45]; 
other exact localized solutions of the BVP for S2 7^ decreasing as x ^ 00 were discussed 
in [46]. 

In the general case, S{x) is not a binomial and the use of RG symmetry (156!) yields an 
approximate analytic solution of the BVP, studied in detail in [44,46] for a Gaussian beam 
with N = exp(— x^). This solution of the BVP allows tracing the evolution of the Gaussian 
beam as the distance from the interface increases, up to a singularity occurring in the 
solution; the singularity has the 2-dimensional structure: both the beam intensity n and 



the derivatives Vx and rix become infinite at the point tf^g'^ = 1/ \/2{a — 13) for a > (3. A 
thorough discussion of this analytic solution and its comparison with the results of other 
approaches (aberration-free approximation [39] and the method of moments [47,48]) were 
carried out in [46]. 



4 Renormalization-group symmetries in nonlocal prob- 
lems of mathematical physics 

The implementation of the RG algorithm in problems of mathematical physics involving 
nonlocal (integral or integrodifferential) relations depends on the form of this nonlocality. 
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On the one hand, the original system of equations can be based on nonlocal relations, as, 
for example, in the kinetic plasma theory, according to which the relation between the 
current density and the charge density in a medium and the distribution function of the 
plasma particles in the Vlasov- Maxwell equations with a self-consistent field is nonlocal. 
The application of the RG algorithm to such nonlocal problems of mathematical physics 
proceeds in accordance with the general scheme in Section [21 the difference is in the 
methods of the calculation of symmetries for nonlocal objects (see [43] and the references 
therein). We note that in the case of problems described by complicated equations, as 
in transfer theory (the Boltzmann integrodifferential equation), only some components 
of the solution or its integral characteristics can have a relatively simple symmetry. For 
instance, in the simplest planar one-velocity transfer problem, the RG invariance is a 
property of the asymptotic form as x ^ oo of the 'density of particles going inside the 
medium' n^{x), which does not feature in the Boltzmann equationjf] 

On the other hand, interesting from the physical standpoint can often be not the 
solution itself over the entire range of its arguments and parameters but some integral 
characteristic, a functional of the solution. For instance, this characteristic can be the re- 
sult of averaging (integrating) with respect to some independent variable or of passing to 
another integral (e.g., Fourier) representation. In this case, we can use the RG algorithm 
to improve the functional of the approximate solution rather than to improve the partic- 
ular solution and the subsequent calculation of the corresponding integral characteristic. 

We now present examples of the implementation of the RG algorithm in problems 
of mathematical physics involving nonlocal relations, which provide illustrations to both 
possible cases. 

4.1 Renormalization-group symmetries for functionals of solu- 



We consider some EVP for local equations and assume that we are interested in an integral 
characteristic of the solution, given by a linear functional of this solution: 



We assume that for a particular solution u of this boundary value problem, the RG 
algorithm has been used to find an RG symmetry with a generator R. Instead of the RG 
transformation group of the solution itself, we are interested in the RG transformation 
group of integral characteristic (1571) . To find an infinitesimal generator of the group, we 
extend the action of the RG symmetry operator R to nonlocal variable For this, we 
represent the operator in the canonical form, that is, make the substitution R ^ Y = xdu-, 
and extend the operator to the nonlocal variable J: 



tions 




(57) 




1 
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The x'^ variable is related to x by means of an integral relation [49] (for brevity, we write 
only one argument of the coordinate of the generator, the one with respect to which the 
integration is performed): 



X 



Considering operator fl59|) in the narrowed space of the variables defining the functional, 
we obtain the required infinitesimal RG symmetry operator for integral characteristic 

dSZD- 



4.1.1 The RG symmetry of functionals in the Hopf equation 

To demonstrate how formulas (!58|) and (1591) actually work for functionals of solutions of 
BVPs, we start with our example of the BVP for the Hopf equation. The algebra of 
RG symmetries of this problem is generated by the three operators in (!23l) . We consider 
the case where we are interested not in the full solution to BVP ( 1251) for all values of 
its arguments and parameters but only in the value at some point of some characteristic, 
which can be defined by a linear functional of form 0571) . For instance, we can be interested 
in the value of the first spatial derivative at x = 0: 

+ 00 

dxu{t, x)\x=o = = ~ J S'{x)u{t, x) . (60) 

— oo 

Using perturbation theory (12T!) in the right-hand side of (!60|) yields the behavior of for 
small t <C 1: 

(uDpt = U! - et [([/°)2 + UiQ)Ul] + O (t^) , [7° ^ d.Ul^^, Ul ^ 9..f/|^^„ . (61) 

To correct the asymptotic behavior of the functional of the solution, which is distorted 
by perturbation theory (1611) . we can use the RG symmetry for (M?]) . As in the derivation 
of solution (125]) . we use the last generator in (12^ in the simplest form, with eip^ = 1, and 
write it in the space of variables {t, e, For instance, for U = x , this operator is 

R = d,-t{u',Yd^o. (62) 

Information about the behavior of the function u^, = 1/ (1 + et) , in the entire range of 
the t variable, including its asymptotic behavior as t oo, can be obtained either with 
the help of (finite) transformations of the group with generator (l62l) . which are similar 
to (^^, or by using the obvious invariant = et — l/u° of RG generator (IB^ . with 
the initial condition M°(t = 0) = 1. We emphasize that we obtain this result without 
explicitly finding solution (125|) . but only using the RG symmetry. Our construction may 
look cumbersome at first glance; however, in more complex problems, the solution is 
typically not known explicitly, but the RG symmetry can be constructed (see, e.g., [50]). 
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4.1.2 The RG symmetry of functionals in quasi-Chaplygin media 

One example of a more complicated situation is the behavior of functionals of solutions 
of the BVP for the equations of quasi-Chaplygin media (135|) and, more specifically, of the 
quantities n{x) and v{x) on the axis x = 0, up to the point where a singularity occurs. 
We claim that this phenomenon can be investigated by applying the RG algorithm to two 
functionals of solutions of BVP (1351) : the density n^{t) = n(t, 0) and the derivative of the 
velocity W'^it) = f^(t,0) calculated on the axis of the beam and related to the solution 
by the formal equalities 

n\t) = J dx6{x)n{t,x) , W°{t) = j dx 6{x)v^{t,x) . (63) 

The boundary conditions for functionals (163|) can be written as 

n°(0) = l, W\0) = 0. (64) 

Although conditions (64) give no information about the dependence of the density n on 
the X variable, such information is incorporated into the RG symmetry operator, whose 
explicit form is determined by the density distribution N{x) for t = 0. We consider 
an example of the problem with the planar geometry, with the 'soliton' profile N{x) = 
cosh~^(x) of the density distribution at the boundary, for the RG symmetry operator as 
in (l38l) . fHOj) . Extending this operator to nonlocal variables (l63l) . we obtain a simpler 
operator in the space {t, n"} [26]: 

R = x""9„o, x"" = 4 - 5n° - + 2(n° - l)r2°n° (n")-^, n° = OtU^ < = dun°. 

(65) 

The RG invariance condition x" = for operator fl65|) leads to an ordinary second-order 
differential equation for the function n^{t), which must be solved with initial conditions 
flMl) and the additional condition {n^/ — 1) j^^Q = 2 for the first derivative, which fol- 
lows from the original equations ( !35l) for a; = 0. This solution t = ^/n^ — reproduces 
the result obtained in (H2l) . but the method is simpler and solution (H2ll is not explicitly 
required. 

We note that the procedure of extending RG generators represented as infinitesimal 
operators of a point group or a Lie-Backhand group is universal, and hence we have a 
common framework for the description of the behavior of characteristics of solutions of 
BVPs ( |T5l) and ( l35l) alike, if we use the reduced description in terms of functionals of 
solutions. 

4.2 Renormalizat ion-group symmetry in the problem of an ex- 
panding plasma bunch 

We now consider the construction of the RG symmetry in the problem where nonlocal 
relations are involved in the definition of the basic manifold. We discuss the problem of an 
expanding plasma bunch in the quasineutral approximation [51]. In this approximation. 
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the dynamics of plasma particles in the planar geometry are determined by solutions of 
kinetic equations for the distribution functions /"(t, x, v) of the different kinds of particles 
(electrons and ions): 

dtr + vdj'' + {eJm^)E{t, x)d,r = (66) 

with additionally imposed nonlocal constraints arising from the conditions of the vanishing 
current and charge density (the quasineutrality conditions): 

/"dt;5^e„r = 0, ldvvJ2^o.r = ^- (67) 

a a 

The electric field is here expressed in terms of moments of the distribution functions: 

E{t, x) = (^J dvv' d,J2 (^/ E ■ (68) 

The initial conditions for system flBBl) . (IB71) correspond to the distribution functions of 
electrons and ions at the instant t = 0: 



n,^, = f^{x,v). (69) 

To construct the RG symmetry, we regard the system of local (166|) and nonlocal (1671) 
equations as the manifold TZAi (step I of the RG algorithm), on which the electric field 
E{t,x) is to be determined. The calculation of the Lie group of point transformations 
admitted by this manifold (step (II) of the RG algorithm) defines a finite-dimensional 
algebra generated by the operators of time and space shifts, the Galilean transformation 
operator, three dilation operators, the quasineutrality operator, and the operator of the 
projective group. The restriction of the group (step (III) of the RG algorithm) on a 
particular solution of problem (jHS]), dSZD, with a spatially symmetric initial distribu- 
tion function with zero mean velocity selects a linear combination of the time translation 
operator and the projective group operator. Under this combination, the approximate 
solution of the initial problem /" = /^(x, f) + 0(t) provided by the perturbation theory 
as t — >■ is invariant, and therefore this linear combination is the RG symmetry operator: 

R = {1 + nH^)dt + nHxd^ + n'^ix - vt)d^ . (70) 

The constant VL can be regarded as the ratio of the characteristic sound speed Cg to the 
initial inhomogeneity scale of the plasma density Lq. 

The invariants of RG operator fITOl) are given by the distribution functions of the par- 
ticles /" and the combinations x / \fl+TV't? and v'^ + Q?{x — vtY . Hence, the construction 
of a solution of the EVP (step (IV) of the RG algorithm) reduces to expressing the dis- 
tribution functions at an arbitrary instant t 7^ in terms of initial data fl69|) with the help 
of these invariants, 

r = /o"(/^"^) , /^"^ = \{^' + ^'(^ - ^tf) + —%{x') , (71) 

/ m„ 
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where the dependence of $o on x' = x/y/l+TWi^ is determined by quasineutrahty con- 
ditions (jnZD- A concrete example illustrating these formulas for a plasma layer formed by 
a group of hot and cold electrons and two kinds of ions can be found in [51]. 

Applications of the RG symmetry operator are not confined to the construction of 
solutions of an initial value problem for Eqns (!66l) . (!67j) or to finding the corresponding 
distribution functions of particles. For practical purposes, more rough characteristics of 
the plasma dynamics is often needed, for instance, the density of the particles (ions) of a 
certain kind n'^{t,x), which can be found by integrating the distribution function: 

oo 

n'^{t,x) = J dvf'^{t,x,v). (72) 

— oo 

Straightforward integration of the distribution function with respect to the velocity cannot 
always be performed analytically because this function may have a complicated depen- 
dence on the invariant 1^°"^ In this case, we can use the extension of the RG symmetry 
operator to a functional of the solution because the density n'^{t,x) is a linear functional 
of f^. The extension of operator (1701) to functional (1721) yields the following operator in 
the narrowed space of the variables {t, x, n'^}: 

R = {i + nH^)dt + nhxd^ - nHnWni . (73) 

The solution of the Lie equations for operator fl73|) with initial conditions fl69l) taken into 
account yields a relation between the invariants of this operator (one of the invariants, 
J3 = x/y/l+TWt^, coincides with the above-mentioned invariant of operator (1701) . and 
the other invariant is J4 = n'^y/T+TW^ at an arbitrary instant t 7^ to their values at 
the initial instant t = 0: J3|t=o= x', Jl\t=o= A/'q(x'). This relation immediately gives 
formulas describing the space-time distribution of the density of the ions of given species 
in terms of their initial density distribution: 

00 

"'^vmm^'ivm^)- - / . (74) 

—00 

We note that the function A/'g also characterizes the energy spectral distribution of ions for 
large times > 1 [51]. Thus, the use of the RG algorithm not only allows constructing a 
solution of problem fl66l) . fl671) . fl69|) for various initial distribution functions of particles [51] 
but also permits finding the law of the evolution of their density and their energy spectrum 
without calculating the distribution functions of the particles explicitly. Similar results are 
obtained not only in the framework of the model of a planar one-dimensional expansion 
but also, for instance, for a spherically symmetric expansion of a plasma bunch [52]. 

5 Conclusion 

We now expound on several important points related to the development and applications 
of the RG algorithm to BVPs of mathematical physics. 
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First of all, we note its universality, meaning that the procedure for the construction 
and the use of RG symmetries is implemented in accordance with the scheme described 
in Section [21 Of course, approaches to the realization of the steps of the algorithm can be 
different depending on the type of the problem under consideration, but the general pat- 
tern of four successive steps remains the same. Our method not merely allows reproducing 
already known solutions in a regular fashion but also produces new solutions. 

Second, the above examples do not exhaust all possible ways of implementing the 
RG algorithm. There is an especially large freedom at the first step, that is, in the 
construction of the original manifold. We have restricted ourself to the description of 
the most typical approaches (extending the list of independent variables, using higher 
symmetries, applying the techniques of approximate symmetries). We left the detailed 
description of the construction of the basic manifold with the use of additional differential 
constraints and methods for the derivation of these constraints on the basis of higher 
symmetries [23] outside the scope of this paper. One special case of a differential relation 
defining a boundary condition is the embedding equation; this is particularly interesting 
in mathematical models based on ordinary differential equations for which the problem 
of symmetry calculation is nontrivial [22-24]. We also left the use of multiparameter 
renormalization groups [22], the construction of approximate RG symmetries involving 
a small parameter in the transformation [44], and integration with respect to the RG- 
transformation parameter [22] without detailed speculation. For a detailed discussion 
of these issues and for applications of RG symmetries, the reader can consult reviews 
[23-26,50,53] and the references therein. 

Third, we note that methods of computer algebra can be used for the construction of 
RG symmetries. In the framework of the general scheme of the RG algorithm, one of the 
central computational procedures is finding a maximal symmetry group of the manifold 
TZAi. Here, it is necessary to derive and solve a system of determining equations, which are 
linear (in the coordinates of the infinitesimal operators) ordinary or partial differential 
equations. This usually amounts to routine calculations, the bulk of which becomes 
quite large for higher-order symmetries and which cannot be performed 'by hand' in a 
reasonably short time; psychologically, this can be a factor constraining the use of the 
RG algorithm. However, by using methods of computer algebra at the second step of the 
algorithm, often allows considerably accelerating the construction of RG symmetries, as 
was shown in the example of the calculation of RG symmetries for equations of quasi- 
Chaplygin media [38]. The prospective gains can be at their greatest if analytic and 
symbolic calculations are combined, when a priori information about the form of the RG 
symmetry extracted from analytic investigations can considerably reduce the time required 
for symbolic calculations. Methods of symbolic calculations can be used for exact and 
approximate RG symmetries alike, which significantly enhances the potentialities of the 
RG algorithm in general. At the same time, analytic approaches used in constructing RG 
symmetries can be helpful in the development of new algorithms for computer algebra 
systems. 

Finally, we indicate possible ways to extend the scope of applications of the RG al- 
gorithm. This can be achieved by covering new objects for which the use of the RG 
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algorithm is not yet standard or by modifying the algorithm itself. 

One example of a new object can be an infinite system of coupled integrodifferential 
equations similar to systems for correlation functions in statistical physics or to systems of 
equations for generalized Green's functions, propagators, and vertex functions in quantum 
field theory. 

As concerns modifications of the algorithm, they are connected in a natural way with 
the general progress in the modern group analysis. This is how it became possible to 
extend the RG algorithm (developed originally for physical problems described by differ- 
ential equations) to nonlocal problems. Certain hopes in this direction are related to the 
progress in group analysis in application to generalized functions [54], further developing 
the theory of approximate symmetries [55], finding new relations between the concept 
of symmetry and conservation laws [56], progress of the theory of partial-invariant solu- 
tions [57,58], and applying group analysis to difference [59] and functional [60] equations. 

This research was carried out with the financial support of the ISTC (grant 2289), 
RFBR (grants 06-02-16103 and 08-01-00291), and LSS-1027.2008.2. 
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